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Currently proposed architectures for long-distance quantum communication rely on networks of 
quantum processors connected by optical communications channels p], Q|. The key resource for 
such networks is the entanglement of matter-based quantum systems with quantum optical fields 
for information transmission. The optical interaction bandwidth of these material systems is a 
tiny fraction of that available for optical communication, and the temporal shape of the quantum 
optical output pulse is often poorly suited for long-distance transmission. Here we demonstrate that 
nonlinear mixing of a quantum light pulse with a spectrally tailored classical field can compress 
the quantum pulse by more than a factor of 100 and flexibly reshape its temporal waveform, while 
preserving all quantum properties, including entanglement. Waveform conversion can be used with 
heralded arrays of quantum light emitters to enable quantum communication at the full data rate of 
optical telecommunications. 

The development of long-distance quantum communication is critical for future quantum cryptography and dis- 
tributed quantum computing applications. Current fiber-optical quantum communication systems rely on direct 
transmission of quantum light pulses, but the attenuation of the fiber imposes a distance limit of tens of kilometers 
for this kind of quantum communication Q by virtue of the no-cloning theorem Q . Quantum repeater architectures 
0, Q promise to circumvent this limit by preparing entangled states over an optical communications channel and 
storing these entangled states as a resource for subsequent quantum communication. Components of quantum 
repeaters have now been demonstrated with a wide variety of physical systems, including single atoms [H-0], atomic 
vapors rare-earth ions in solids [ill, [HI, quantum dots (Hf, and NV-centres [3]. The common thread among 

these demonstrations is the manipulation of quantum light pulses by matter-based quantum emitters. The temporal 
waveform of such emitters is typically a single-sided exponential with decay constant on the order of 1 nanosecond, 
which cannot be readily mode-matched to the smooth, broadband pulses desirable for telecommunications, posing 
a substantial disadvantage for quantum-emitter approaches to quantum networking. Recent attempts to overcome 
these issues include increasing the emitter bandwidth to the GHz range by the use of nonresonant interactions 
[il| . shaping the temporal waveform by placing the emitter in a nonlinear resonator [l6j], temporal modulation of 
single- photon wavepackets fl7l . Il8j , and nonlinear frequency conversion experiments with single photons [l9l - |22| . 

Here we present an efficient and straightforward method of quantum optical pulse shaping and compression that im- 
mensely simplifies the interfacing of quantum emitters with telecommunications networks. Our method preserves the 
full quantum statistics of the input field, including entanglement and any other multimode correlations, while enabling 
compression by more than a factor of 100, along with flexible reshaping of the temporal waveform. In particular, our 
method enables time/ wavelength transduction of spontaneously emitted photons from quantum emitters into short, 
smooth pulses at telecommunications wavelengths. As shown in Fig. [TJ the input field undergoes three-wave mixing 
(3WM) with a frequency-chirped classical laser pulse. For an appropriate choice of classical laser intensity and chirp, 
the 3WM product radiation has the same spectrum as the desired target mode, but receives the quantum statistics 
of the input. The 3WM output is then dechirped with a second pulse shaper to match the temporal wavefunction of 
the target mode. The waveform converter extends the classical time-lens technique, which has achieved remarkable 
results in compressing and stretching classical pulses [H,[23|, to the quantum domain and to arbitrary pulse reshaping. 

We describe the 3WM process using slowly varying bosonic field operators 1 i>j(z,t) with j = 1 the input mode and 
j = 2 the mode generated by 3WM. Here z measures the distance along the propagation axis in a frame comoving at 
the group velocity, and t measures the duration of the interaction between the three fields. (This coordinate convention 
is to be contrasted with the classical nonlinear optics convention in which z measures the interaction length in the 
3WM medium and t is the time of arrival at the detector.) The escort laser pulse contains > 10 10 photons and can 
be approximated as a classical field that remains unaffected by 3WM. For an escort pulse much longer than the other 



escort laser 




FIG. 1: Schematic of the quantum optical waveform converter. A nonclassical light input (originating, e.g., from a quantum 
emitter in a high-finesse optical resonator) is combined with a highly chirped classical pulse. The combined fields undergo 
three-wave mixing (3WM) in a nonlinear crystal, transferring the spectral modulation of the classical pulse onto the 3WM 
output. The output is separated from the original fields and passed through a pulse shaper to remove the residual phase, 
producing a pulse in the target mode that inherits the quantum state of the input mode. The colours under the pulse envelope 
represent the frequency variation during the pulse length, with the variation of colours greatly exaggerated for clarity. 

pulses, the Hamiltonian becomes (see Methods) 

H 3WM = ihfl J dz e l ^ z) ^\t> 2 + h.c. (1) 

where f2 is determined by the nonlinear coupling constant and the intensity of the escort pulse, and where cj>(z) is the 
phase of the escort field. The quantum field operators evolve as 

*i(z,i) = cosfii*i(z,0) + sinme^ (z) * 2 (z,0) (2) 
*2(*,*) = -sinOte^ (z) *i(z,0) + cosfi** 2 (*,0) (3) 

If the fields leave the 3WM medium after an interaction time T = 7r/(20), the solution for mode 2 at times t > T is 
just ^2(2, T) = — e~ l ^ z ^ \{z, 0). The quantum state of mode 1 is perfectly transferred into mode 2, while mode 2 
has acquired the phase —<j){z) + tt from mode 3. To temporally match the 3WM output to the target pulse shape, 
the output pulse shaper then removes the undesired relative phases of the spectral components, performing a unitary 
transformation on the output field operator. The spontaneous emission from a quantum emitter, with a single-sided 
exponential waveform, can be converted into a much shorter Gaussian pulse by choosing (see Methods) 

0( z ) = [ d( eirUe-^^h (4) 
a Jo 

where r is the spontaneous emission lifetime and the target amplitude is proportional to e~ z /( 2 ( ctr ) ). The phase 
modulation of Eq. (d} is visualised in Fig. [5J 

The ideal quantum waveform conversion described above will not be achieved with unit fidelity in real 3WM media 
because of dispersion. We now show that the fidelity F can nevertheless exceed 99.9% for readily achievable experi- 
mental parameters. For pure input states, F — |(V>ideal I ^disp)| ; where the result of ideal evolution is written jV'ideai) 
and the result with dispersion included is IV'disp)- If the input system is entangled with another quantum system, 
the fidelity of the final entangled state is simply the average fidelity of the eigenstates of the input density operator, 
weighted by their corresponding eigenvalues. A perturbative analysis of the dispersive evolution (see Methods) shows 
that the error is dominated by mismatch between the group velocities t>i,2,3 of the input, output, and escort fields, 
which we parametrise by v = (vi — «2)/2 and v e = V3 — (vi + v%)/2. For a pure input wavefunction A(z) with 
characteristic length scale L, we define dimensionless velocities u — v/vq, u e = v e /vo, where vq = QL/(2ir). The error 
can be minimised by adding the compensation phase 



A op t(z) = -(u e - u)4>'(z)L 



(5) 
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FIG. 2: a) Escort phase modulation function <j>(z) for conversion from a single-sided exponential waveform to a Gaussian 
waveform with compression ratio t/ct = 100. b) Visualisation of the spectral chirp dcj>{z)/dz imposed on the initial waveform 
by 3WM. Height indicates waveform amplitude, colour indicates local frequency after escort phase imprinting. The hue of the 
colour is proportional to d<j>(z)/dz as calculated from Eq. @. 

to the initial escort phase 4>(z). For an average photon occupation (n), we then obtain 



The ratio u e /u is set by the crystal dispersion alone, but u varies with the escort laser intensity I csc as u oc I esc . 
Thus the fidelity can be made arbitrarily close to 1 by increasing the escort laser power. Eq. ((6]) can then be rewritten 
as 1 — F op t = (u/u crT ) 2 for some it crr < 1. It can be seen from Eq. (j6]) that the perturbation theory breaks down for 
a pulse with an arbitrarily sharp leading edge; a perturbative analysis in momentum space shows that 1 — F opt oc u 
in this limit. In practice, the time required to excite a quantum emitter is never exactly zero, so the leading edge of 
the pulse is smoothed over the excitation timescale. We also perform a full numerical simulation for a single-mode 
single-photon input state (see Methods). This confirms that the effect of group- velocity dispersion is insignificant in 
all cases of interest. 

Figure [3] shows the analytic and numerical error estimates for two cases of particular experimental interest. Case 1: 
the conversion of 370 nm photons from a Yb + ion Q to the 1550 nm telecommunications band using periodically 
poled lithium niobate, for which u e /u « —2/3. The simulated error closely follows the perturbative result for small 
values of u up to a compression ratio t / a — 200. For an escort laser pulse of energy ~ 1/iJ and duration of 150 ns 
(> 20 ryb+) in a 50 mm long crystal waveguide, one finds u = 0.013 and error of 1 — F = 7 x 10~ 4 at compression 
ratio of 100. Case 2: the conversion of 780 nm photons from a Rb atom to the telecommunications band, for which 
one can arrange u e /u = —1 by poling lithium niobate for type II phase matching (25|. Here the error is much lower 
for similar escort laser parameters, v < 10 5 m/s for any choice of output wavelengths in the telecommunications 
band, so u < 10~ 4 and 1 - F « 10~ 4 . 

We have shown that three-wave mixing with a modulated classical field can reshape and compress the waveform of 
a quantum light pulse while faithfully maintaining the quantum information carried by the photons. A quantum 
light pulse produced by a quantum emitter with a lifetime of nanoseconds can be converted to a Gaussian pulse 
with a duration of tens of picoseconds that is compatible with standard telecommunications protocols. The low 
error of waveform conversion is compatible with schemes for fault-tolerant quantum communication over long 
distances [l[. Quantum waveform conversion enables simultaneous time- and wavelength-division multiplexing of 
the pulses from an array of quantum emitters up to the limit of channel capacity, massively increasing quantum 
communications bandwidth. Current DWDM systems with 50 GHz channel spacing achieve their maximum capacity 
for transform-limited pulses with ~ 20 ps duration, while the dispersive effects of long-haul fibre transmission require 
the pulses to have a smooth temporal waveform. As each pulse arrives from the emitter array, it can be simultaneously 
converted to this ideal waveform and sorted into an appropriate DWDM channels. The rate of entangled pair genera- 
tion in a quantum network is then limited only by the telecommunications bandwidth and the size of the emitter array. 




(6) 



Methods 



Quantum 3WM Hamiltonian We analyze the 3WM process using the canonical quantization method [26-28] . For simplicity, 
the electric field polarisation vectors are assumed to lie along the 3WM crystal axes, as for a periodically poled or otherwise 
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FIG. 3: Error induced by quantum waveform conversion from a single-sided exponential pulse of time constant r and rise time 
0.02r to a Gaussian pulse of 1/e 2 time constant a. Errors are less than 1CF 3 for readily achievable experimental parameters 
(see text), a) Error at compression ratio r/a — 100 as a function of dimensionless group- velocity mismatch u for conversion 
of 370 nm photons to 1550 nm in lithium niobate. Solid line: perturbative prediction. Points: simulation results. Dashed 
line: Best fit of 1 — F = (u/u clr ) 2 to simulation points, b) Error scale u CTT as a function of compression ratio r/cr. Solid line: 
perturbative prediction. Points: simulated values computed from least-squares fits to simulation results. As expected, the 
simulation results match well to the perturbative theory, c), d) are the same as a), b), but for conversion of 780 nm photons to 
1550 nm in type-II matched lithium niobate, near the special point v — —v e at which the escort phase term in Eq. <(6j vanishes. 
Errors are even lower than for conversion of 370 nm photons, but at compression ratio above 10 the perturbation theory breaks 
down and higher-order GVM dominates the error. 



noncritically phase-matched crystal. We assume conservation of momentum and energy for the carrier waves and retain only 
phase-matched processes. The x' 2 ' nonlinear Hamiltonian is then 

H 3WM = iha / dz *1*2*3 + h.c. (7) 



where the coupling constant a is determined by the material nonlinear susceptibility and the beam geometry and is readily calcu- 
lated in the classical limit. The escort field is taken to be a classical field of constant intensity that is phase-modulated to impart 
the desired spectral modulation to mode 2. We write ^3 = £exp[i0(z + v e t)], where £ is the (real, positive) amplitude of the 
classical escort field and </>(z + v e t) is the local phase. With the definition SI = a£, Eq. ([7]) reduces to Eq. {TJ for the case v e — 0. 

Phase functions for waveform shaping 

For simplicity, we assume that the input mode and the desired target mode are both transform-limited. To match the power 
spectrum of the 3WM product to the desired spectrum, 4>( z ) should satisfy |cb(fc)| = 1/ dz ct\ (z)e II *' z - ) e _l ' cz | , where cb(fc) is 

the Fourier transform of 02(2). In general, one can satisfy this constraint by numerical least-squares minimization. However, 
when the input and target bandwidths differ substantially, the method of stationary phase applies to the integral and 4>{z) has 
a closed-form solution in this limit when the target is Gaussian. Writing 02 (k) tx e~ k ^ 2rj \ we find 



a 



Ci 

a + b I d( 2 a 2 (( 2 ] 



(8) 



where erf is the error function, f~ 1 (x) = y denotes the solution of f(y) = x, and the constants a,b are set by the boundary 
conditions of the transformation. After 3WM, the phase of ai (k) is nontrivial and the 3WM product pulse is therefore not 
transform-limited. The output pulse shaper applies a spectral compensation phase 7(fc), implementing the unitary transforma- 
tion 

*o„,(z) = ^ J dz<t 2 (z y t = T) J dke^ k) e^ ( - z) (9) 
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Choosing 



7(fe) = -<f> 



removes the unwanted phase, so that the output pulse is transform-limited with the desired spectrum. When a\(z) oc e 



(10) 



V(2m 2 ) 



is also Gaussian, equations ((HJ and (|10|) reduce to • 
the solution is that given in Eq. ((4]). 



(z) = az /(2/i) and j(k) = —fife /(2a), while for a single-sided exponential 



Dispersive evolution and error in state transfer 

Errors in the state transfer arise from group-velocity mismatch between the three modes in the 3WM medium. In the comoving 
frame with velocity v = (vi +V2)/2, the quantum Hamiltonian for group- velocity mismatch (GVM) and group- velocity dispersion 
(GVD) can be written as [29| 



H, 



disp 



= E 

3=1,2 



dz 



ihvi 



— 3 -Vi 



2 d 



+ 



9*, 



4 dz dz 



+ h.c. 



(11) 



while the 3WM Hamiltonian |TJ is also modified because the phase 4>{z) of the escort field phase evolves under dispersion (see 
Supplementary Discussion). Moving to the interaction picture with respect to the original 3WM Hamiltonian ((TJ), one derives 
the additional unitary evolution due to dispersive effects, Uaisp- A second-order Dyson series solution for Udi Bp shows that 
GVM mixes the vacuum noise of the initially unoccupied mode 2 into the state transfer, while GVD has a negligible effect. 
The removal of phase by the output pulse shaper just implements a unitary transformation on the 3WM output field, which 
has no effect on the fidelity. The fidelity is then evaluated as F = ({7dis P )|, where the expectation value is taken with respect 
to the initial states of modes 1 and 2 and any other systems entangled with mode 1. Compensating the phase according to Eq. 
© is found to minimise the error independently of the input state. Taking an initial pure state in mode 1 and the vacuum 
state in mode 2, we obtain Eq. ((6}. 

Numerical simulations of error in state transfer 

The Heisenberg equations of motion for the field operators 9 2 are linear, so the operator of the target field after time 
t will be a linear combination of the initial field operators. If the input pulse has a single spatial mode A(z), such that 

|V(0)> = /[o t ] |0) = / [/ efcc4(z)*J(z)] |0) = J2 n c « I 71 )- the quantum state at time t is 



M0)> = / 



/^{^(M)^(,) + A 2 (M)*1(,)} 



(12) 



with Ai(z,0) — A(z) and Ai{z,0) = 0. The A n obey the same linear equations as \l/ ni but are c-number amplitudes rather 
than operators. Simulating these equations allows us to calculate the fidelity as 



F = 



]T| Cn | 2 j dzA*(z)e- l4 ' {z) A 2 (z,T) 



(13) 



For fidelities close to unity and a correctly compensated phase, this reproduces the linear dependence on (n) found in the 
perturbative calculation ©. For definiteness, we only show results for a single- photon input state. 
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Supplementary discussion: Perturbative calculation of waveform conversion error 



Dispersive evolution 



Error in the state transfer arises from dispersion in the 3WM evolution. To calculate the error, we treat all dispersion- 
related terms as perturbations to the nondispersive 3WM Hamiltonian 



H 3WM = ihn J dz^\(z,t)^ 2 (z,t)e im +h.< 



(14) 



which induces ideal state transfer by the unitary evolution Uq(T) over the time T = 7r/(20). Mode 2 is assumed to 
be initially unoccupied throughout the calculation, as for the applications discussed in the main paper. 

We consider dispersion up to second order, i.e., group velocity mismatch (GVM) between modes 1, 2, and 3, and 
group velocity dispersion (GVD) within each of these modes. In nondegenerate mixing, as considered here, GVM has 
much larger effects than GVD. The GVD and GVM parameters are expressed using the variables v = (v% — v 2 )/2 and 
v e = v 3 — (vi + v 2 )/2 for GVM (where Wj = duji/dki) and /3, = d 2 uii/dk 2 for GVD. We define e ^ 1 the perturbative 
expansion parameter and take v,v e ~ 0(e), Pi ~ 0(e 2 ). The dispersion of the quantum modes 1 and 2 is described 
by the Hamiltonian 



1 >^i{z,t) |-^^ 2 ( z ,t) + h.c. 



dz 



dz 



dz 

d 2 «f\{z,t) 

dz 2 



d 2 n(z,t) 
dz 2 



V 2 (z,t)+h.c. 



(15) 
(16) 



The dispersion of the escort laser (mode 3) affects the quantum state transfer indirectly through the 3WM Hamiltonian. 
At 0(e 2 ), the escort mode is governed by the evolution equation 



d^ 3 
dt 



df 3 ijkf^s 

'■ dz 2 dz 2 



(17) 



As will be seen, the error can be minimised by precompensation of the phase function. For an escort pulse of constant 
intensity, "J^ oc e~ lx ^ z,t \ we write the compensated phase function at t — as x(z, 0) = 4>{z) + A(z). The solution of 
Eq. (HU is 



X(Z, t) = <t>( Z - V e t) + A(Z - V e t) - ?f [(<j)'(z)) 2 - i<j>"{z)] 



(18) 



where the prime applied to functions (as in (f>') indicates differentiation with respect to z. Including the effect of 
dispersion, the 3WM Hamiltonian becomes 



^3WM = #3WM 



ihQ, I dz^ 



{(z,£)V 2 (z,t) (e*< 



We go to an interaction frame with respect to 7?3wm to obtain the perturbation Hamiltonian 

Vj(t) = U^ (H disp + H^ M - H 3WM )U « Vx{t) + V 2 (t) 



(19) 



(20) 



where V%(t) — 0(e) and V 2 (t) = 0(e 2 ). The Dyson series gives the unitary operator Ui(t) describing the perturbed 
evolution in the interaction frame: 



Uj(t) = i-l J* dt Vj(t) - 1 J* df J* dt" ViWViW) 

Expanding the Dyson series to 0(e 2 ) gives 



Un{t) 



U cff (t) a 1 + U n (t) + U n (t) + U 12 (t) 
ft 



(21) 



(22) 



dtVx{t) U 2 i(t) 



dtV 2 (t) U 12 (t) 



dt' / dt" Vi{t')Vx(t") (23) 
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It can be shown that V\ and V 2 are both Hermitian, so the expectation values of Uu(T) and U2i(T) are purely imagi- 
nary. Un(T) contributes to the error only as 0(v 2 ,v 2 ) ~ C(e 2 ), showing that we must include U\2(T) in a consistent 
perturbative expansion of the error to lowest order in GVM. Moreover, U^iiT) contributes only as 0(f3f) <~ C(e 4 ). 
We consider only error terms up to C(e 2 ), so we discard U2i(T). Hence GVD has no effect on the error in this analysis. 

The solution of the nondispersive Hamiltonian lets us express ^i(z, t) and ^ 2 (z,t) in terms of $(z) = ^\(z,t = 0) 
and T(z) = ^ 2 (z,t = 0), giving 



dz 



(p. sin 2flt A + (v sin 2 fit - v e Qt sin 2Qt) $ f $ — iv cos 2flt — — $ 
-2fie** cos 2^ (A - v e t<t>') + <e**t; sin 2Clt f^T - $ f 



+ h.c. 



(24) 



Since Vi(£) is normally ordered and mode 2 is initially unoccupied, terms that involve only mode 2 do not affect the 
fidelity and arc omitted from the expressions. We wish to evaluate the error at the end of the state transfer, which 
occurs at time t = T = 7r/(2fi), giving 



Un{T) = ~i} I M*)* f (*)*(*) 

h{z) = 4ftA(z) + n(v - v e )<f>'(z) 



(25) 
(26) 



so that A ~ 0(e) or higher if A is to minimise the error in state transfer. The remaining term in U c s(T) is computed 
to be 



U 12 (T) 



1 



dz 1 dz 2 {~2h{z 1 )h(z 2 )^{z 1 )^{z 1 )^{z 2 )^{z 2 ) 

-8« 2 C0S [ ( j>(z 1 )-<f>(Z2)}^(z 1 )<f>'(Z2)^(z 1 )r(z 1 )r^Z2Mz2) 



_2 Ue i (0(zi)-</Hz2)) 



9+(2l) , t(2l) x (2 , ) (^l 4(S2) -xt (a) ^) 



*^T(*0-*n*i)^W)*(*) 



-8w 2 e i(0(zi)_0(z2)) 



v y T(zi)T T (z 2 ) v + $ T (zi) , v / $(z 2 ) 



cfei dz 2 



d,Z2 



dz\ 



dz 2 



(27) 



Here g±(z) = 47rftA(z) — (7r 2 ± 8)v e (j>'(z) and we have eliminated terms with purely imaginary expectation values, as 
these terms do not contribute to the error at 0(e 2 ). 



Fidelity calculation 

We quantify the error in the state transfer by computing the fidelity F between the actual output state of the waveform 
converter and the ideal dispersion-free output state. The final pulse shaping just implements a unitary transformation 
on the 3WM output field, which has no effect on the fidelity For many applications one wishes to convert entangled 
states involving both mode 1 and some other modes, so the full state before waveform conversion takes the form 



|*<> = E 



*j \4 j) ) 



\x U) ) 



(28) 



where the ) are orthonormal states of mode 1 and the |x^) are orthonormal states over the other systems. 
Writing Uf u \\(T) as the unitary evolution under the full dispersive Hamiltonian, we have 



0) 



U) 



(29) 
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Thus, since we can always find the transfer fidelity for entangled input states by taking an appropriate weighted sum 
over pure-state fidelities, we need only calculate the transfer fidelity of a pure input state \ipi}- In many cases, such a 
state is characterised by a mode creation operator 

at = / dzA{z)&{z) 



(30) 

where $(z) — ^>i(z,t = 0) and the mode wavefunction A(z) is normalised as J dz \A(z)\ 2 — 1. An initial fc-photon 
number state in the mode is given by (a^) k |0) and a general single-mode initial pure state — X^^Lo ° k 1^) can 
be written as \tpi) = f(a')\0), where |fc) denotes a number state and f(x) is defined through the series expansion 

The pure-state fidelity is given by the squared overlap between the state obtained from phase-compensated dispersive 
evolution and that obtained from uncompensated nondispersive evolution. In an interaction frame with respect to 
the nondispersive 3WM Hamiltonian, we have 



F = 



4>i 



Ul s (T) 



4>i 



(31) 



and using the matrix elements derived in the Appendix, we find 



(iPi\U n (T)\^i 
(4>% I Ui 2 (T) | ipi 



tin) 

i 

64ft 2 



dz h(z)\A(z)\' 



(32) 



d Zl dz 2 {-2{n(n - l))h{z 1 )h{z 2 )\A{z 1 )\ 2 \A{z 2 )\ 2 + 2(n)h(z 1 )h(z 2 )A*(z 1 )A(z 2 )S(z 1 - z 2 

8(n)v 2 e cos [0( Zl ) - <f>{z 2 )] ^{z 1 )4> l {z 2 )A*{z 1 )A(z 2 )8{z 1 - z 2 ) 
i{n)ve^M-*^ [g+izt) (A* (z 1 )A(z 2 )d Z2 S(z 1 - z 2 ) - A*(z 1 )A'(z 2 )S(z 1 - z 2 )) 

-g-(z2) (A'*(z 1 )A(z 2 )5(z 1 - z 2 ) - A*(z 1 )A(z 2 )d Zl S(z 1 - z 2 ))] 
8{n)v 2 e i ^~^ [A'*(z 1 )A'(z 2 )S(z 1 - z 2 ) + A* (z 1 )A(z 2 )d Zl d Z2 S(z 1 - z 2 ) 

-A*(z 1 )A'(z 2 )d zl 6(z 1 - z 2 ) - A* (zi)A(z 2 )d Z2 8{zi - z 2 )]} 



64ft 2 



2(n(n-l))( I dzh\A{ 2 ) +2(n) I dz h 2 \A\ 2 - 8(n)v 2 I dz\(j>' A\ 2 



+8(n)vv e I dz (j)' [A*(iA + <p'A) + h.c] - 8(n)v 2 I dz \2A' - i(j>'A\ 2 



^2 - !)) ^ dz h \ A \^j + («) J dz h2 \ A \ 2 - 4 (") J dz \ 2vA ' - { ( V + Ve)<P'A\ 



(33) 

(34) 
(35) 



where n is the number operator of the mode, (n) denotes the expectation value of n, and we have retained only real 
terms. The fidelity is optimised when h(z) — so that 



Aopt(s) 



n(v e - v)(f>'{z) 
4ft 



Fopt = 1 ~ 8ft^ ' dz \ 2vA '^ ~ ^ v + v e)4>'{z)A{z)\ 



(36) 
(37) 



which is just Eq. ([5]). 
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Appendix: Matrix elements for pure states 

To evaluate the matrix elements involved in the fidelity calculation, we first observe that 

[$(*), at] = Jd(A(0 [$(*), $t( C )] = J d( A(()6(( - z) = A(z) 

[a,a^} = jdQA*(C) [*(C),a + ] = Jd(\A(()\ 2 = l 

[$(«), (a*)*] = fcA(z)(a t ) fc " 1 
(0|a J '(a t ) fe |0) = fc!5 jfe 

Writing (n) for the expectation value of the photon number in we find 

oo oo ^. 

Mi = I ^1)^2) |V») = EE ^7m <° I a^i)^)^ I 0) 



fc=0 .;=() 



= A*(z 2 )A(zi) £ ^fc 2 (0 | a fe " V)*" 1 | 0) 



fc=0 



fc! 



(38) 

(39) 

(40) 
(41) 

(42) 
(43) 



= A*( Zl )A(z 2 ) k\c k \ 2 (44) 

fc=0 

= (n)A*( Zl )A(z 2 ) (45) 

OO OO ^, 

M 2 = | & {zi )$(z 2 )&(z 3 )$(z 4 ) |^) = EE ^ft (° I o^i)*^)*^)*^)^)* | 0) (46) 

fe =oj=o 

= A*(^M( Z4 ) E E -3g=jk (0 I a'- 1 ^)*^)^)*- 1 | o) (47) 
fc =0j =o v^'^- 

= A*( Zl )A(z 4 ) E E "TsM* (° I a '" 1 (* f (*3)*te) + <J(*S - z 3 )) (at)^ 1 I 0) (48) 
k=o 3 =o v k -3- 

00 I 1 2 

= A*( Zl )A(z 4 ) -^-k 2 ((fc - l) 2 (fc - 2)\A*(z 3 )A(z 2 ) + 6( Z2 - z 3 )(fc - 1)!) (49) 



fc=0 



(n(n - l))A*( Zl )A(^)A*( Z3 )A( Z4 ) + (n)A*( Zl )A(z 4 )S(z 2 - z 3 ) 



Then 



$t( zl ) 



d$(z 2 ) 



$ t (z 1 )$(z 2 )$ t (z 3 ) 



dz 2 

rf$t( Zl ) rf$( Z2 ) 
dzi d Z2 
d<&( Zi ) 

d Z4 



(50) 

(51) 
(52) 

^ = 9 Z4 M 2 = (n(n - l))A*(z 1 )A(z 2 )A*(z 3 ),4'(z 4 ) - (n)A*(z 1 )A'(z 4 ),5(z 2 - z 3 ) 

(53) 



^=d tl Mi = (n)A*( Zl )A'(z 2 ) 
4>i) = d Zl d Z2 M! = (n)(A*Y( Zl )A'(z 2 ) 



and similarly for the other matrix elements involving mode 1. For mode 2, which is initially unoccupied, we calculate 

(54) 



M 3 = (0\ T(z!)Tt(z 2 ) | 0) - 5( Zl - z 2 ) 
dT( Zl ) 



, -Tt(z 2 ) 
dz\ 

dX{z Y )dT\z 2 ) 
dz\ dz 2 



= d zl M 3 = d zl S( Zl - z 2 ) 
) = d Zl d Z2 5{z - z 2 ) 



(55) 
(56) 



